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ABSTRACT 

We explore the infrared spectrum of a three-dimensional dynamical model of planet HD 209458b as a function 
of orbital phase. The dynamical model predicts day-side atmospheric pressure-temperature profiles that are much 
more isothermal at pressures less than 1 bar than one-dimensional radiative-convective models have found. The 
resulting day-side thermal spectra are very similar to a blackbody, and only weak water absorption features are 
seen at short wavelengths. The dayside emission is consequently in significantly better agreement with ground- 
based and space-based secondary eclipse data than any previous models, which predict strong flux peaks and deep 
absorption features. At other orbital phases, absorption due to carbon monoxide and methane is also predicted. 
We compute the spectra under two treatments of atmospheric chemistry: one uses the predictions of equilibrium 
chemistry, and the other uses non-equilibrium chemistry, which ties the timescales of methane and carbon monox- 
ide chemistry to dynamical timescales. As a function of orbital phase, we predict planet-to-star flux ratios for 
standard infrared bands and all Spitzer Space Telescope bands. In Spitzer bands, we predict 2-fold to 15-fold 
variation in planetary flux as a function of orbital phase with equilibrium chemistry, and 2-fold to 4-fold variation 
with non-equilibrium chemistry. Variation is generally more pronounced in bands from 3-10 than at longer 
wavelengths. The orbital phase of maximum thermal emission in infrared bands is 15-45 orbital degrees before 
the time of secondary eclipse. Changes in flux as a function of orbital phase for HD 209458b should be observable 
with Spitzer, given the previously acheived observational error bars. 

Subject headings: planetary systems, radiative transfer, binaries:eclipsing, stars: individual (HD 209458) 



1. INTRODUCTION 

Astronomers and planetary scientists are just beginning to 
understand the atmospheres of the short period giant planets 
known as "hot Jupiters" or "Pegasi planets." A key subset of 
these planets are those that transit the disk of their parent stars, 
which make them well-suited for follow-up studies. The most 
well studied of these transiting hot Jupiters is the first to be dis- 
covered, HD 209458b (Charbonneau et al. 2000; Henry et al. 
2000), a 0.69 Mj planet that orbits its Sun-like parent star at a 
distance of 0.045 AU. 

Currently, considerable work on hot Jupiters is occuring on 
both the observational and theoretical fronts. In the past few 
years, several groups have computed dynamical atmosphere 
models for HD 209458b in an effort to understand the struc- 
ture, winds, and temperature contrasts of the planet's atmo- 
sphere (Showman & Guillot 2002; Cho et al. 2003; Cooper 
& Showman 2005; Burkert et al. 2005; Cooper & Showman 
2006). If the planet has been tidally de-spun and become locked 
to its parent star, dynamical models are surely needed to under- 
stand the extent to which absorbed stellar energy is transported 
onto the planet's permanent night side. With the launch of the 
Spitzer Space Telescope, there is now a platform that is well- 
suited for observations of the thermal emission from hot Jupiter 
planets. 

Spitzer observations spanning the time of the planet's sec- 
ondary eclipse (when the planet passes behind its parent star) 
have been published for HD 209458b at 24 pm (Deming et al. 
2005), TrES-1 at 4.5 and 8.0 pm (Charbonneau et al. 2005), 



and HD 189733b at 16^m (Deming et al. 2006a). In all cases, 
the observed quantity is the planet-to-star flux ratio in Spitzer 
InfaRed Array Camera (IRAC), InfraRed Spectrograph (IRS), 
or Multiband Imaging Photometer for Spitzer (MIPS) bands. 
The dual TrES-1 observations are especially interesting because 
they allow for a determination of the planet's mid-infrared spec- 
tral slope. Several efforts have also been made from the ground 
to observe the secondary eclipse of HD 209458b. Although 
no detections have been made, some important, occasionally 
overlooked upper limits at K and L band have been obtained. 
These include Richardson et al. (2003a) around 2.3 /im, Snellen 
(2005) in K band, and Deming et al. (2006b) at L band. These 
ground-based observations constrain the flux emitted by the 
planet in spectral bands where water vapor opacity is expected 
to be minimal; therefore, emitted flux should be high. 

This influx of data has spurred a new generation of radiative- 
convective equilibrium models, whose resulting infrared spec- 
tra can be compared with data (Fortney et al. 2005; Burrows 
et al. 2005; Seager et al. 2005; Barman et al. 2005; Fortney 
et al. 2006; Burrows et al. 2006). See Marley et al. (2006) 
and Charbonneau et al. (2006) for reviews. The majority of 
these models are one-dimensional. Authors weight the inci- 
dent stellar flux by 1/4, to simulate planet-wide average condi- 
tions, or by 1/2, to simulate day-side average conditions (with 
a cold night side). Barman et al. (2005) have investigated two- 
dimensional models with axial symmetry around the planet's 
substellar-antistellar axis and computed infrared spectra as a 
function of orbital phase. Iro et al. (2005) have extended one- 
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dimensional models by adding heat transport due to a simple 
parametrization of winds to generate longitude-dependent tem- 
perature maps, but they did not compute disk averaged spec- 
tra for these models. Very recently Burrows et al. (2006) have 
also investigated spectra and light curves of planets with var- 
ious day-night effective temperature differences, assuming ID 
profiles for each hemisphere. These one- and two-dimensional 
radiative-convective equilibrium models have had some success 
in matching Spitzer observations, but Seager et al. (2005) and 
Deming et al. (2006b) have shown that ground-based data for 
HD 209458b do not indicate promiment flux peaks at 2.3 and 
3.8 /im, which solar composition models predict. 

The various dynamical models for HD 209458b (Showman 
& Guillot 2002; Cho et al. 2003; Cooper & Showman 2005; 
Burkert et al. 2005; Cooper & Showman 2006) are quite var- 
ied in their treatment of the planet's atmosphere. We will not 
review them here, as that is not the focus of this paper. In 
general, temperature contrasts in the visible atmosphere are ex- 
pected to be somewhere between 300 and 1000 K. What has 
been somewhat lacking for these dynamical models are clear 
observational signatures, which would in principle be testable 
with Spitzer or other telescopes. The purpose of this paper is 
to remedy that situation. Here we generate infrared spectra and 
light curves as a function of orbital phase for the Cooper & 
Showman (2006) (hereafter: CS06) dynamical simulation. We 
present the first spectra generated for three-dimensional models 
of the atmosphere of HD 209458b. 

2. METHODS 

Here we take the first step towards understanding the effects 
of atmospheric dynamics on the infrared spectra of hot Jupiters. 
A consistent treatment of coupled atmospheric dynamics, non- 
equilibrium chemistry, and radiative transfer would be a con- 
siderable task. In a coupled scheme, given a three-dimensional 
P-T grid at a given time step, with corresponding chemical 
mixing ratios, the radiative transfer scheme would solve for 
the upward and downward fluxes in each layer. These fluxes 
would be wavelength dependent and would differ from layer to 
layer. The thermodynamical heating/cooling rate, which is the 
vertical divergence of the net flux, would then be calculated. 
The dynamics scheme would then use this heating rate, together 
with the velocities and P-T profiles in the grid at the previous 
time step, to a calculate the chemical abundances, velocities, 
and P-T profiles on the grid at the new time step. The pro- 
cess steps forward in time, as the radiative transfer solver again 
finds the new heating/cooling rates. The emergent spectrum of 
the planet could be found at any stage. In our work presented 
here, we performed a simplified calculation that contains some 
aspects of what will eventually be included in a fully consistent 
treatment. 

2.1. Dynamical Model 

Our input pressure-temperature (P-T) map is from CS06. As 
the dynamical simulations are described in depth in Cooper & 
Showman (2005) and CS06, we will only give an overview 
here. The CS06 model employs the ARIES/GEOS Dynam- 
ical Core, version 2 (AGDC2; Suarez & Takacs 1995). The 
AGDC2 solves the primitive equations of dynamical meteorol- 
ogy, which are the foundation of numerous climate and numer- 
ical weather prediction models (Holton 1992; Kalnay 2002). 
The primitive equations simplify the Navier-Stokes equations 
of fluid mechanics by assuming hydrostatic balance of each ver- 



tical column of atmosphere. Forcing is due to incident flux from 
the parent star, through a Newtonian radiative process described 
in CS06. The CS06 model is forced from the one-dimensional 
radiative-convective equilibrium profile of Iro et al. (2005), 
which assumes globally averaged planetary conditions. 

For their simulations, CS06 take the top layer of the model 
to be 1 mbar. The model atmosphere spans ~15 pressure scale 
heights between the input top layer and the bottom boundary at 
3 kbar. CS06 use 40 layers evenly spaced in log pressure. A P— 
T profile is generated at locations evenly spaced in longitude (in 
5 degree increments) and latitude (in 4 degree increments). The 
72 longitude and 44 latitude points create 3 168 P-T profiles. In 
Figure 1 we show the CS06 grid at 3 pressure levels, near <~2, 
20, and 200 mbar. Previous work has shown these levels likely 
bracket the pressures of interest for forming mid-infrared spec- 
tra (Fortney et al. 2005). All three panels of Figure 1 use the 
same brightness scale for easy comparison between pressure 
levels. 

At the 2 mbar level, where radiative time constants are short 
(Iro et al. 2005; Seager et al. 2005), the atmosphere responds 
quickly to incident radiation. Winds, though reaching a speed 
of up to 8 km s" 1 , are not fast enough to lead to significant devi- 
ations from a static atmosphere, which implies that the hottest 
regions remain at the substellar point. The arrows indicate a 
wind pattern that attempts to carry energy radially away from 
the hot spot. At this pressure, the atmosphere somewhat resem- 
bles the two-dimensional radially symmetric static atmosphere 
of Barman et al. (2005), who found a hot spot at the substellar 
point and a uniformly decreasing temperature as radial distance 
from this point increased. The night side appears nearly uni- 
form, and colder. 

At the 25 mbar level, it is clear that a west-to-east circulation 
pattern has emerged at the equator, and the center of the planet's 
warm region has been blown downstream by ^35 degrees. The 
wind from the west dominates over the predominantly radially 
outward wind seen at the 2 mbar level. Day-night temperature 
contrasts are not as large at this pressure as they are at 2 mbar. 

At the 220 mbar level, the center of the hot spot has been 
blown downstream by ^60 degrees. This jet extends from the 
equator to the mid-latitudes; the gas in the jet is warmer than 
gas to the north or south. The radiative time constants become 
longer the deeper one goes into the atmosphere (Iro et al. 2005). 
Hence, winds are better able to redistribute energy, leading to 
weaker temperature contrasts, which cannot simply be charac- 
terized as "day-night." 

2.2. Logistics & Radiative Transfer 

Unlike other published models, we stress that the spectra 
generated here are from a dynamical atmosphere model that is 
not in radiative-convective equilibrium. Each of the 3168 P-T 
profiles from CS06 have, without modification (aside from in- 
terpolation onto a different pressure grid), been run through our 
radiative transfer solver. No iteration is done to achieve radia- 
tive equilibrium. The equation of radiative transfer is solved 
with the two-stream source function technique described in 
Toon et al. (1989). This is the same infrared radiative trans- 
fer scheme used in Fortney et al. (2005), Fortney et al. (2006) 
and M. S. Marley et al. (in prep.). We ignore contributions due 
to scattered stellar photons, as discussed below. 

At a given orbital position, the CS06 map in longitude is re- 
mapped into an apparent longitude (as seen from Earth), while 
the latitude remains unchanged. See Figure 2 for a diagram. 
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Here we ignore the 3.4 degrees that the orbit differs from being 
exactly edge-on (Brown et al. 2001). At a given orbital angle 
f — for each patch of the planet — the cosine of the angle 6 from 
the sub-observer point, /j,, is calculated. This \x is consistently 
included when solving the radiative transfer, which means that 
the effects of limb darkening (or brightening) are automati- 
cally incorporated. We interpolate in a pressure-temperature- 
abundance grid (Lodders & Fegley 2002), such that any given 
point in the three-dimensional model has the atomic and molec- 
ular abundances that are consistent with that point's pressure 
and temperature. 

At any given time, the one-half of the planet that is not visible 
is not included in the radiative transfer. The 1584 visible points, 
at which the emergent specific intensity (erg s" 1 cm" 2 Hz" 1 sr" 1 ) 
is calculated, are then weighted by the apparent visible area of 
their respective patches. These intensities are summed up to 
give the total emergent flux density (erg s" 1 cm" 2 Hz" 1 ) from 
the planet. The spectra generated by the patch-by-patch version 
of the code were tested against spectra from one-dimensional 
gray atmospheres and our previously published hot Jupiter pro- 
files. 

Emergent spectra are calculated from 0.26 to 325 /xm, but 
since we ignore the contribution due to scattered stellar flux, 
the spectra at the very shortest wavelengths have little mean- 
ing. However, for the radiative-equilibrium HD 209458b model 
published in Fortney et al. (2005), they found that scattered stel- 
lar flux is greater than thermal emission only at wavelengths 
less than ~0.68 /im. We expect that a considerable amount of 
"visible" light that may eventually be seen from hot Jupiters 
is due to thermal emission. We note that by 1 \im thermal 
emission is 100 times greater than scattered flux. Here we will 
present spectra for wavelengths from 1 to 30 /zm. 

We note that the radiative transfer at every point is solved in 
the plane-parallel approximation. While this treatment is suf- 
ficient for our purposes, we wish to point out two drawbacks. 
The first is that, near the limb of the planet, we will tend to 
overestimate the path lengths of photons emerging from the at- 
mosphere, as the curvature of the atmosphere is neglected. The 
second issue also occurs near the limb. We cannot treat pho- 
tons whose path, in a completely correct treatment, would start 
in one column but emerge from an adjoining column. How- 
ever, atmospheric properties in any two adjoining columns are 
in general quite similar. The former issue, that of the plane- 
parallel approximation, is likely more important, and should be 
addressed at a later time, when data precision warrants it. Here 
we note that in our tests 95% of planetary flux emerges from 
within 75 degrees of the sub-observer point, such that these 
limb effects will have little effect on the disk-summed spectra 
and light curves that we present. 

When generating our model spectra, we use the elemental 
abundance data of Lodders (2003) and chemical equilibrium 
compositions computed with the CONDOR code, as described 
in Fegley & Lodders (1994), Lodders & Fegley (2002), and 
Lodders (2002). In §2.3 we will discuss deviations from equi- 
librium chemistry in the atmosphere of HD 209458b, as calcu- 
lated by CS06. At this time, we ignore photochemistry, which 
has been shown by Visscher & Fegley (2005a) to be reasonable 
at P > 10 mbar. For the most part, the infrared spectra of hot 
Jupiters are sensitive to opacity at 10 < P < 200 mbar (Fortney 
et al. 2005), so equilibrium chemistry calculations are probably 
sufficient. As discussed in Fortney et al. (2006), we maintain 
a large and constantly updated opacity database, which is de- 



scribed in detail in R. S. Freedman & K. Lodders (in prep.). The 
pressure, temperature, composition, and wavelength-dependent 
opacity is tabulated beforehand using the correlated-k method 
(Goody et al. 1989) in 196 wavelength interval bins. The re- 
sulting spectra are therefore of low resolution. However, low 
resolution is suitable for the task at hand, as we are inter- 
ested in band-integrated fluxes and the radiative transfer must 
be solved at 1584 locations on the planet at many (here, 36) 
orbital phases. 

In Williams et al. (2006), which focused on examining asym- 
metrical secondary eclipse light curves caused by dynamical re- 
distribution of flux, we investigated the effects of limb darken- 
ing for HD 209458b for the CS06 map. These effects are hard 
to disentangle from general brightness variations due to tem- 
perature differences generated by dynamics. Limb darkening 
in a particular wavelength would be manifested as a brightness 
temperature that decreases towards the limb relative to a bright- 
ness temperature map computed assuming normal incidence at 
every point. On the planet's day side, where the P-T profiles 
are somewhat isothermal, limb darkening is not expected to be 
significant. Indeed, only at angles greater than ^80 degrees 
from the subsolar point was day-side limb darkening as large as 
100-200 K calculated. As noted, due to our plane-parallel ap- 
proximation, this is likely to be somewhat of an overestimation. 

2.3. Clouds and Chemistry 

Before calculating spectra for this dynamic atmosphere, it 
is worthwhile to step back and look at the atmospheric P- 
T profiles with an eye towards understanding the effects that 
clouds and chemistry may have on the emergent spectra. In 
Figure 3, we plot a random sampling of the 3168 P-T profiles 
and compare them to cloud condensation and chemical equi- 
librium boundaries. It should be noted that there is a greater 
density of profiles on the left side of the plot. These are profiles 
from the relatively uniform and cool night hemisphere. Also 
of note is that at pressures less than ^200 mbar, the warmer 
profiles are fairly isothermal, with quite shallow temperature 
gradients. 

The CS06 profiles are nearly everywhere cooler than the con- 
densation curves of iron and Mg-silicates. These clouds will 
form, but for these profiles, cloud bases would lie deep in the 
atmosphere at pressures greater than 1 kbar, as this is the high- 
est pressure at which the profiles cross the condensation curves 
of these species. The opacity from such deep clouds would have 
no effect on the spectrum of HD 209458b. 

The temperatures of the CS06 simulation are computed as 
departures from the equilibrium temperature profile of Iro et al. 
(2005), as described in CS06. If a warmer base profile had 
been selected, such as Fortney et al. (2005) or Barman et al. 
(2005) (which are 100-300 K warmer at our pressures of in- 
terest), these curves could be shifted to the right by 100-300 
degrees. Marley et al. (2006) provide a graphical comparison 
of profiles computed by Iro et al. (2005), Fortney et al. (2005), 
and Barman et al. (2005) under similar assumptions concern- 
ing redistribution of stellar flux and atmospheric abundances; 
they find differences of up to 300 K at 100 mbar. These differ- 
ences can probably be attributed to different opacity databases, 
molecular abundances, radiative transfer methods, and perhaps 
incident stellar fluxes. For the Fortney et al. (2005) profile, the 
major heating species are H2O bands from 1 to 3 /im and neu- 
tral atomic Na and K, which absorb strongly in the optical. The 
major cooling species are CO, at 5 /im, and H2O in bands at 
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3 and 4-10 /im, although for a colder profile, such as CS06 
night side profiles, CH4 would also cool across these wave- 
lengths. CS06 chose the Iro et al. (2005) profile because Iro et 
al. also computed atmospheric radiative time constants, which 
other authors have not done to date. Here, the computed night 
side profiles cross the condensation curve of Na2S at low pres- 
sure. Showman & Guillot (2002) and Iro et al. (2005) pointed 
out that the transit detection of weak neutral atomic Na absorp- 
tion by Charbonneau et al. (2002) could be explained if a large 
fraction of Na was tied up in this condensate. 

If the predictions of equilibrium chemistry hold for this at- 
mosphere, then it would appear that at pressures up to 100 
mbar, CH4 would be the main carbon carrier on the planet's 
night-side, while CO would be the main carrier on the day 
side. This would lead to dramatically different spectra for 
these hemispheres, especially at near-infrared and mid-infrared 
wavelengths. 

However, equilibrium chemical abundances would be ex- 
pected only if chemical timescales are faster than any mixing 
timescales. The relative abundances of CO and CH4 can be 
driven out of equilibrium if mixing timescales due to dynamical 
winds are faster than the timescale for relaxation to equilibrium 
of the (net) reaction 

CO + 3H 2 ^CH 4 + H 2 0. (1) 

The time constant for chemical relaxation toward equilibrium 
is a strong function of temperature and pressure; it is short 
in the deep interior and extremely long in the observable at- 
mosphere. Deviations from CH4/CO equilibrium are due to 
the long timescale for conversion of CO to CH4, and are well 
known in the atmospheres of the Jovian planets (Prinn & Bar- 
shay 1977; Fegley & Prinn 1985; Fegley & Lodders 1994; 
Bezard et al. 2002; Visscher & Fegley 2005b) and brown dwarfs 
(Fegley & Lodders 1996; Noll et al. 1997; Griffith & Yelle 
1999; Saumon et al. 2000, 2003, 2006), where vertical mixing 
can be due to convection and/or eddy diffusion. CS06 follow 
timescales for CH4/CO chemistry taken from Yung et al. (1988) 
and find that vertical winds of 5-10 m sec -1 push this system 
out of equilibrium at pressures of interest for the formation of 
infrared spectra. They find that the quench level, where the 
mixing timescale and chemistry timescale are equal, is ~ 1-10 
bar, which is below the visible atmosphere. Above this quench 
level, the mole fractions of CO and CH4 are essentially con- 
stant. They find that the CH4/CO ratio becomes homogenized 
at pressures less than 1 bar everywhere on the planet. For their 
nomical models this CH4/CO ratio is ^0.014, meaning that the 
majority of carbon is indeed in CO. However, a non-negligible 
fraction of carbon remains in CH4 on both the day and the night 
hemispheres. We note that CS06 ignore possible effects due 
to photochemistry on the CO and CH4 mixing ratios. Given 
the few explorations into hot Jupiter carbon and oxygen photo- 
chemistry to date (Liang et al. 2003, 2004; Visscher & Fegley 
2005a), it is unclear how important photochemistry will be in 
determining the abundances of these species at pressures of tens 
to hundreds of millibar. CS06 did explore other effects such 
at atmospheric metallicity and temperature. If the atmosphere 
is greater than [M/H]=0.0, or if the atmosphere is hotter, the 
CH4/CO ratio would be even smaller. See CS06 for additional 
discussion. 

In our spectral calculations, we find that significant differ- 
ences arise depending on our treatment of chemistry. We will 
therefore investigate the effects of a few chemistry cases, as 
explained below, and shown in Table 1. We label our equi- 



librium chemistry trial "Case 0," as it is the standard case. 
For "Case 1," we fix the CH4/CO ratio at 0.014 (as found by 
CS06) at all temperatures and pressures along the profiles. Con- 
sistently incorporating the increasing CH4/CO ratio at depth 
(P > 1 bar) would be difficult with previously tabulated opaci- 
ties, and would have little to no effect on the emergent spectra. 
Recall that 23% of available oxygen is lost to the formation of 
Mg-silicate clouds (Lodders 2003), which in this model have 
cloud bases near 1 kbar. The remaining oxygen is almost en- 
tirely found in CO and H2O. As the CH4 and CO abundances 
are fixed, for Case 1 we will also fix the abundance of H2O; 
the mixing ratios of CH4, CO and H2O are consistent with the 
amount of available oxygen at T~ 1200 K and pressures of tens 
of millibars. We also briefly consider a "Case 2," in which the 
CH4/CO ratio has been further reduced, in an ad-hoc manner, 
which could be due to quenching of the abundances at a hot- 
ter temperature or the photochemical destruction of CH4. The 
CH4/CO ratio is reduced by nearly a factor of 500, to 3 x 10~ 5 . 
The mixing ratios of CO and H2O are nearly the same as in 
Case 1, although the slightly increased CO abundance (due to 
the drop in the abundance of CH4 and conservation of carbon 
atoms) uses up some oxygen at the expense of H2O. This case is 
valuable for comparison purposes because it shows the spectral 
effects of a negligible CH4 abundance. In all cases, the mix- 
ing ratios of all other chemical species are given by equilibrium 
values. However, only CH4, CO, and H2O have a discernable 
impact on the spectra. 

In Figure 4 we plot mixing ratios of CH4, CO, and H2O pre- 
dicted from equilibrium chemistry along three P-T profiles in 
the atmosphere of HD 209458b. In gray are the mixing ratios 
for the one-dimensional profile of Fortney et al. (2005). The 
abundance of CH4 is negligible. In thin black and thick black 
are the mixing ratios at the sub-stellar point and the anti-stellar 
point, respectively, of the CS06 simulation. The dominant car- 
bon carrier is clearly CO, rather than CH4, except at P <100 
mbar for the anti-stellar point. Most of the nightside has a 
chemistry profile similar to the anti-stellar point, and hence, 
strong CH4 absorption will be seen. At the top of the plot, ar- 
rows indicate the assigned mixing ratios of Cases 1 and 2 for 
these three molecules. 

In §3 we make quantitative comparisons to current ground- 
based and space-based infrared data. While we find that the 
agreement is good for the model presented here, we wish to 
stress that this study is exploratory. This is the first study that 
quantitatively explores the influence of atmospheric dynamics 
on the emergent spectra of a hot Jupiter atmosphere. The great- 
est uncertainty likely lies in the calculation of heating/cooling 
rates with a simple Newtonian cooling scheme in the dynam- 
ical model, as discussed in Cooper & Showman (2005) and 
CS06. Another issue is that CS06 calculate temperature devia- 
tions from a one dimensional radiative-convective profile pub- 
lished in Iro et al. (2005). However, our emergent spectrum is 
calculated using our radiative transfer solver, which is a differ- 
ent code, and we use different abundances (Lodders 2003) than 
these authors (Anders & Grevesse 1989). However, we find that 
when solving the radiative transfer for the Iro et al. profile we 
obtain a T e a that differs by only 1 % from their value. 

In their non-equilibrium chemistry work, CS06 chose the 
abundances of Lodders & Fegley (1998), while here we use 
(Lodders 2003), in order for the most natural comparison with 
our previous work (Fortney et al. 2005, 2006). Obviously no 
choice would be fully self consistent. However, the choice of el- 
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emental abundances is certainly a smaller concern than the cur- 
rent debate concerning the correct chemical timescale for con- 
version of CO to CH4 in planetary atmospheres (see Yung et al. 
1988; Fegley & Lodders 1994; Griffith & Yelle 1999; Bezard 
et al. 2002; Visscher & Fegley 2005b). In addion, CS06 have 
shown that a 300 K decrease in temperature, for instance, leads 
to a 20-fold increase in CH 4 /CO (see §5.3), so temperature un- 
certainties likely swamp any abundance issues. We stress that 
our focus in the following sections is on highlighting the spec- 
tral differences between a model that accounts for dynamical 
redistribution of energy around the planet and one that does 
not. Indeed the magnitude of the spectral effect suggests that 
additional, more internally self consistent work, is clearly ap- 
propriate. 

3. INFRARED SPECTRA 
3.1. Spectra as a Function of Orbital Phase 

We now turn to the predicted infrared spectrum of our dy- 
namical model atmosphere. As has been shown by many au- 
thors since Seager & Sasselov (1998), the infrared spectra of 
hot Jupiters are believed to be carved predominantly by absorp- 
tion by H2O, CO, and, if temperatures are cool enough, CH4. 
In general, absorption features of hot Jupiters are predicted to 
be shallower than brown dwarfs of similar effective tempera- 
ture and abundances. This is predicated on hot Jupiters hav- 
ing a significantly shallower atmospheric temperature gradient, 
which is due to the intense external irradiation by the parent 
star. Figure 5 shows absorption cross sections per molecule for 
CH4, CO, and H2O from 1 to 30 /urn. To avoid clutter on our 
spectral plots, we will not label absorption features, so referring 
to Figure 5 will be helpful. 

Our first spectral calculation is for Case chemistry, as a 
function of orbital phase. This is shown in Figure 6. We show 
the emitted spectrum at six orbital locations: zero (during the 
transit), 60, 120, 180 (during the secondary eclipse), 240, and 
300 degrees, where the degrees are orbital degrees after transit. 
At zero degrees (the black spectrum), we see the full night side 
of the planet. The temperature gradient is fairly steep, leading 
to deep absorption features. Most carbon is in the form of CH4, 
leading to strong methane absorption in bands centered on 1 .6, 
2.2, 3.3, and 7.8 fim. Absorption due to H2O at many wave- 
lengths is also strong. Absorption due to CO at 2.3 and 4.5 
is quite weak, due to its small mixing ratio. 

Sixty degrees later (the red spectrum), the hemisphere we 
see is 2/3 night and 1/3 day. The part of the day hemisphere 
we see has the hot west-to-east jet coming towards us. This 
will later be an interesting point of comparison with the 300 
orbital degree spectrum (in magenta), which is also 2/3 night. 
Here all absorption features are muted because, due to the jet 
shown in Figure 1, more than ~l/3 of the planet is showing the 
nearly isothermal profiles representative of the day side. These 
profiles lead to a blackbody-like spectrum. At 120 orbital de- 
grees (in green), we are seeing 2/3 day side, and that part of the 
day side is the hottest, due to the wind from the west blowing 
the hot spot downstream. The absorption features have become 
very weak. 

The peak in total infrared flux actually occurs 153 orbital de- 
grees after transit. During the secondary eclipse, 27 degrees 
later, the planet is less luminous. During the secondary eclipse 
(in dark blue) the planetary spectrum is essentially featureless 
long-ward of 1.8 /jm. The fully illuminated hemisphere pos- 
sesses profiles that are nearly isothermal in our pressure range 



of interest. Intriguingly, this leads to a spectrum that is very 
similar to a blackbody. For comparison, we plot a dashed black 
curve that is the thermal emission of a 1330 K blackbody. It 
plots behind the 180 degree spectrum for all A > 4 /zm. 

As we move to orbital phases after the secondary eclipse, we 
see that the spectra are not symmetric with their reflections on 
the other half of the orbit. The light blue spectrum is signifi- 
cantly lower than the green one, even though at both points we 
are seeing 2/3 day and 1/3 night. This is again due to the west- 
to-east jet. At 240 orbital degrees, the hottest part of the day 
side is not in view. At 300 orbital degrees (in magenta) 1/3 of 
our view is the coolest part of the day side, along with 2/3 of 
the night side. The spectrum is quite similar to what was seen 
during the transit. The planet appears least luminous 3 1 orbital 
degrees before transit. 

3.2. Spectra with Alternate Chemistry 

Since CS06 find significant deviations from equilibrium 
carbon chemistry, it is useful to examine spectra with non- 
equilibrium chemistry. We now investigate the spectra using 
Cases 1 and 2. Figure 7 shows comparison spectra at two or- 
bital phases: zero and 180 orbital degrees, which is during tran- 
sit and secondary eclipse, respectively. On the day side, we see 
that the emergent spectrum is essentially exactly the same, in- 
dependent of chemistry. Because the P—T profiles are nearly 
isothermal, essentially no information about the chemistry can 
be gleaned from the spectrum. For comparison with our pre- 
vious work, in Figure 7, we plot with a dotted line the cloud- 
less planet-wide average spectrum of HD 209458b from Fort- 
ney et al. (2005), which shows absorption due to H2O and CO 
and strong flux peaks between these absorption features. 

We will now examine the small inset in Figure 7. This is 
a zoom-in of the secondary-eclipse spectra near 2.2 ^m. The 
three thick horizontal bars connected by dashes make up an in- 
teresting observational one-sigma upper limit that was first re- 
ported in Richardson et al. (2003a). These data were analyzed 
again and presented in Seager et al. (2005). The measurement 
was a relative one; only the vertical distance between the cen- 
tral band and two side bands is important: they can be moved 
up or down as a group. Specifically, the flux in the central band 
cannot exceed the combined flux in the two adjacent bands 
by an amount greater than vertical distance shown (which is 
0.45 x 10~ 6 erg s" 1 cm" 2 Hz" 1 ). Our dynamical day-side spec- 
tra, which are nearly featureless due to the nearly isothermal 
atmosphere, easily meet this constraint. We note that an up- 
per limit obtained by Snellen (2005) in K band, which covers 
a similar wavelength range, has an error bar that is too large to 
distinguish among published models. 

In contrast to the models presented here that include dynam- 
ics, published ID models to date all predict large flux differ- 
ences between 2.2 and the surrounding spectral regions 
when abundances are solar. The ID models with solar abun- 
dances presented in Seager et al. (2005), for example, were 
not able to meet the 2.2-^m flux-difference constraint, because 
the H2O bands adjoining the 2.2-/im flux peak were too deep. 
Seager et al. (2005) suggested that if C/O > 1, then little H 2 
would exist in the atmosphere, and the observational constraint 
could be met. Here, we propose instead that dynamics produces 
a relatively isothermal dayside temperature, in which case the 
constraint can easily be met even with solar C/O. 

Nevertheless, the observational constraint is not yet firm 
enough to fully rule out the ID radiative-equilibrium models 
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at solar abundances. The Fortney et al. (2005) solar-abundance 
models, for example, marginally meet the current observational 
flux-difference constraint. (Fortney et al. 2005 find a flux dif- 
ference of 0.43 x 10" 6 erg s" 1 cm" 2 Hz" 1 .) The cause for the 
discrepancy between the ID models in Fortney et al. (2005) 
and Seager et al. (2005) remains unclear but could occur if the 
P-T profiles in Seager et al. were steeper than those in Fortney 
et al., leading to deeper bands in the former study. More pre- 
cise flux-difference observations like that of Richardson et al. 
(2003a) during secondary eclipse would help confirm or rule 
out the ID solar-abundance models. Even a factor of two de- 
crease in the observational upper limit between the "in" and 
"out" bands could rule out the solar-abundance Fortney et al. 
(2005) model, while supporting the isothermal-dayside models 
(e.g., CS06) and the C/O > 1 models with weak water bands. 
Alternatively, the detection of a small flux peak at 2.2 fim would 
lend support to the ID models with C/O < 1, while constraining 
the temperature gradient on the planet's dayside. 

The fact that both the solar-abundance CS06 models (pre- 
sented here) and the radiative-equilibrium ID models with C/O 
> 1 produce almost no flux difference between 2.2 /im and 
the surrounding continuum implies that the 2.2-(im band can- 
not be used to distinguish between these alternatives. Instead, 
additional constraints at other wavelengths will be necessary to 
discriminate between them. Measurements of flux differences 
surrounding CO (rather than H2O) bands could provide such 
a test. The radiative-equilibrium C/O > 1 models, although 
lacking water bands, would presumably have strong CO bands 
and would hence predict a strong flux difference between the 
center of a CO band and the surrounding continuum. On the 
other hand, because the dayside atmosphere is nearly isother- 
mal in the CS06 models, these circulation-altered models would 
predict little flux difference between the CO band and the sur- 
rounding regions even in the presence of large quantities of CO. 
An observation of minimal flux difference across CO bands as 
well as H2O bands would support isothermal-dayside models 
like CS06 while arguing against the radiative-equilibrium ID 
models with steep temperature gradients. 

Although observations of TrES-1 may not be directly ap- 
plicable to HD 209458b, as TrES-1 may have a r eff 300 K 
cooler, there is an important issue to note for our discussion 
here. The models of Fortney et al. (2005), Seager et al. (2005), 
and Barman et al. (2005) show a mid-infrared spectral slope 
that is bluer than observed from 4.5 to 8.0 /um by Charbonneau 
et al. (2005) for TrES-1, although the model of Burrows et al. 
(2005) appears to be consistent with this slope. Showman & 
Cooper (2006) and Fortney et al. (2006) have discussed that 
models with an atmospheric temperature inversion would give 
infrared spectra with a redder spectral slope, due to molecular 
emission features, in better agreement with observations. The 
CS06 model of the day side of HD 209458b, which lacks the 
strong negative temperature gradient of these radiative equilib- 
rium models, also leads to a redder mid-infrared spectral slope. 

On the planet's night side, we see significant spectral differ- 
ences between the three chemistry cases. In Cases 1 and 2, the 
CH4 mixing ratio is constrained to a small abundance, weak- 
ening these absorption features. Absorption due to CH4, CO, 
and H 2 is readily seen in Case 1. Case 2 shows essentially 
the same CO and H2O absorption, but CH4 absorption is no 
longer seen. The transition, as a function of orbital phase from 
deep to essentially nonexistent absorption features in Cases 1 
and 2, are similar to what was seen in Case 0. In the interest of 



conciseness, and because Case 2 is ad-hoc, hereafter we con- 
centrate on Cases and 1 . The spectra for Cases 1 and 2 are 
essentially the same, except in regions of CH4 absorption; we 
will highlight these differences when necessary. It is important 
to remember that in Case 1 and Case 2 chemistry, the mixing 
ratios of our principal absorbers, CH4, CO, and H2O are fixed. 
Therefore, changes in the spectra with orbital phase are only 
due to changes in the planetary P-T profiles on the visible disk, 
due to the rotation of the planet through its orbit. 

One can integrate the spectrum of the planet's visible hemi- 
sphere over all wavelengths, as a function of orbital phase, to 
determine the apparent luminosity of the planet at all phases. 
Here we divide out 4irR 2 a to calculate the apparent effective 
temperature (T e ff) of the visible hemisphere. This is plotted in 
Figure 8 for Cases and 1 . We can clearly see that in both 
cases, the time of maximum flux precedes the time of secondary 
eclipse by ^27 degrees, or 6.3 hours. Since the spectra of the 
two cases overlap around the time of secondary eclipse, so do 
the plots of T e ff. At other orbital phases, the Case curve al- 
ways plots lower. The largest effect is before the time of transit. 
The effect is tied to the CH4 abundance. If CH4 is able to at- 
tain a large mixing ratio, it leads to an atmosphere that has a 
higher opacity, meaning one can not see as deeply into the at- 
mosphere. One then reaches an optical depth of <~1 higher in 
the atmosphere, which is significantly cooler for night-side P-T 
profiles, leading to a lower T e ff. The point of minimum plane- 
tary flux precedes the transit by 3 1 to 37 degrees, depending on 
the chemistry. 

In their (very similar) previous dynamical model, Cooper & 
Showman (2005) predicted a time of maximum planetary flux 
of 60 degrees (or 14 hours) before secondary eclipse. The large 
timing difference between that work and this one is due al- 
most entirely to the choice of "photospheric pressure" made in 
Cooper & Showman (2005). They chose a pressure of 220 mbar 
in their work, which is deeper than the "mean" photosphere 
that we find here. At higher pressure, the radiative timescales 
are longer, such that winds are able to blow the atmosphere's 
hottest point farther downstream. Cooper & Showman (2005) 
and Showman & Cooper (2006) previously discussed how their 
prediction varied as a function of the chosen photospheric pres- 
sure. 

4. INFRARED LIGHT CURVES 
4.1. Spitzer Bands 

For the foreseeable future, the most precise data for under- 
standing the atmospheres of hot Jupiters will come from the 
Spitzer Space Telescope. For HD 209458b, only an observation 
at 24 /jm (the shortest wavelength MIPS band), has been pub- 
lished. It seems likely that observations in all four IRAC bands 
(3.6, 4.5, 5.8, and 8.0 /im), as well as IRS at 16 ^m, will be ob- 
tained within the next year or two. As such, we have integrated 
the spectra of our planet models and a Kurucz (1993) model 
of star HD 209458 over the Spitzer bands in order to generate 
planet-to-star flux ratios as a function of orbital phase. These 
are plotted for Cases and 1 in Figure 9. The stellar model fits 
the stellar parameters derived in Brown et al. (2001) and is the 
same model used in Fortney et al. (2005). 

The behavior of the planet-to-star flux ratios is quite interest- 
ing. While the T e ff of the planet was found to reach a maximum 
27 degrees before secondary eclipse, the behavior in individual 
bands is more varied. For instance, in both chemistry cases, 
the planetary flux in the 24 fim band peaks only 15 degrees be- 
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fore secondary eclipse. This is because the "photospheric pres- 
sure" is at a lower pressure in this band than the planet's "mean 
photospheric pressure." As previously discussed, at lower pres- 
sures, the radiative time constants are shorter, and the atmo- 
sphere is able to more quickly adjust to changes in incident flux. 
At higher pressure, winds are better able to blow the planet's 
hot spot downstream. One should keep in mind that the light 
curves generated are a function of the dynamical calculation 
and the radiative transfer. It is the radiative transfer calculation 
that determines how deeply into the atmosphere (and therefore, 
to what temperature) we see. 

The 3.6 /mi band peaks earliest, —27 degrees before transit, 
as the T e ff does as well. This band shows a 15-fold variation 
in flux (peak to trough) as a function of orbital phase for Case 
because it encompasses a significant CH4 absorption feature 
that waxes and wanes. Since the abundances of CH4 and CO 
are not free to vary in Case 1, the flux ratios in this case do not 
show the large amplitudes found in some bands in Case 0. The 
dotted lines in the Case 1 panels are for Case 2 in the 3.6 and 
8.0 /im bands, where CH4 absorption occurs. The flux varia- 
tion in these bands is even further reduced as CH4 absorption 
is not seen. (See Figure 7.) At 24 /im, our model is 1.6 sigma 
higher than the secondary eclipse data point published by Dem- 
ing et al. (2005), indicating that the planet may be dimmer at 
180 orbital degrees than we predict. In all bands, differences 
in chemistry between the two cases have essentially no effect 
on the timing of the maxima in planetary flux. However, since 
the night side is much more sensitive to chemistry, the minima 
can vary by as much as 20 degrees in bands that are sensitive to 
CH 4 absorption. 

4.2. Standard Infrared Bands 

The results of ground-based observations of flux from hot 
Jupiters have been mixed. All searches for visible light have 
yielded only upper limits (Charbonneau et al. 1999; Collier 
Cameron et al. 2002; Leigh et al. 2003), which have ruled out 
some models with bright reflecting clouds. In the near infrared, 
specifically for HD 209458b, Richardson et al. (2003b,a) have 
constrained molecular bands of CH4 and H2O. The constraint 
on emission at 2.3 /im between H2O absorption features was 
shown in Figure 7. The predicted planet-to-star flux ratio really 
does not become favorable until wavelengths longer than 3 /im, 
which may continue to challenge observers. 

In Figure 10, we plot planet-to-star flux ratios in the Mauna 
Kea Observatory (MKO) H (~1.6 /mi), K (~2.2 /mi), L' (-3.8 
/mi), and M' (~4.7 /mi) bands. Since the wavelength ranges 
of the L' and M' bands have significant overlap with the IRAC 
3.6 and 4.5 /mi bands, the predicted ratios are quite similar. In 
dark blue, we show the V band upper limit of -0.0007±0.0014 
from Deming et al. (2006b). The model is just above the 1 
sigma error bar. This is a better match than is attained with 
one dimensional radiative equilibrium models (Deming et al. 
2006b), which predict a flux peak just short of 4 /im, as shown 
in Figure 7. The Fortney et al. (2005) HD 209458b radiative 
equilibrium model, which has somewhat muted flux peaks com- 
pared with similar models by other authors, has a planet-to-star 
flux ratio of 0.00114 in L' band. In dotted blue in the Case 1 
panels is our L' band (which encompasses CH4 absorption) pre- 
diction for Case 2 chemistry. Looking at shorter wavelengths, 
in the H and K bands, the planet is predicted to be dimmer, 

5 As described in Deming et al. (2006b), this observation was actually performed 
planet-to-star flux ratio at secondary eclipse increases by 4% when using this nam 
are unchanged. 



while the star is brighter, leading to low flux ratios. The peak 
emission in these bands does occur earlier than in the Spitzer 
bands. For instance, the H band peak is 42 orbital degrees be- 
fore secondary eclipse. 

5. DISCUSSION 
5.1. Effective Temperature 

For the CS06 dynamical model of the atmosphere of HD 
209458b, we find that the apparent T e s of the visible hemisphere 
is strongly variable, with a maximum of 1 390 K and a minimum 
of 915 K for equilibrium chemistry and 1025 K for (the proba- 
bly more realistic case of) disequilibrium chemistry. This leads 
to a luminosity of the visible planetary hemisphere that varies 
by factors of 5.3 and 3.4, respectively for these two cases. For 
the one dimensional Iro et al. (2005) planet-wide profile, we 
derive a T e s = 1325 K, 1% less than found by the authors. The 
day-side T e ff for the dynamical model is not significantly larger 
than this planet-wide T e ff. As an additional comparison, we cal- 
culate the mean of the planetary luminosity over the entire or- 
bit and then convert to T e ff to find a mean T e ff for Cases and 
1. This gives T e ffS of 1195 K and 1227 K, respectively, show- 
ing that the CS06 model, at the pressures levels that radiate 
to space, is as a whole somewhat colder than the one dimen- 
sial profile of Iro et al. (2005). This is likely a consequence of 
the radiative forcing scheme employed in CS06, which will be 
reinvestigated when models that consistently couple radiative 
transfer and dynamics are developed. 

The ~400 K T e ff contrasts found are predominantly due to 
changes in the temperature structure of the hemisphere that is 
visible as function of orbital phase. In addition, atmospheric 
opacity makes an important contribution when the CH4/CO ra- 
tio is free to vary with pressure and temperature. When the 
CH4 mixing ratio is below that predicted by equilibrium chem- 
istry, this leads to a lower opacity atmosphere, for a given P-T 
profile. While the few individual bands of CO are somewhat 
stronger than those of CH4, CH4 absorption across the planet's 
broad the 2-10 /im flux peak dominates over the two bands of 
CO. One is able to see more deeply into a CH4-depleted atmo- 
sphere, leading to a higher T^. 

5.2. HD 209458b Infrared Data 

Currently, the only published secondary eclipse datum from 
Spitzer for HD 209458b is the 24 /im detection of Deming et al. 
(2005). The models presented here have a planet-to-star flux 
ratio during secondary eclipse that is 1 .6 sigma larger than this 
observational data point. Together with our excellent agreement 
with the ground-based data at 2.3 /im from Richardson et al. 
(2003a) and Seager et al. (2005), and our 1.1 sigma difference 
with the 3.8 /im data from Deming et al. (2006b), we regard 
this as excellent agreement — significantly better than has been 
previously obtained with one dimensional radiative equilibrium 
models. 

It is interesting to discuss a few issues that arise if the flux 
ratios are indeed 10-25% less than we calculate here, as indi- 
cated by the L' and 24 /im band data. For instance, perhaps 
day-night temperature contrasts in the atmosphere are not as 
large as predicted by CS06, leading to smaller deviations from 
a "planet- wide" T e ff — 1200K. This may involve radiative time 

a narrow band centered on 3.8 £tm, within the standard L' band. Our calculated 
band, compared to standard L'. Since this is a small correction, our conclusions 
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constants that are significantly longer than predicted by Iro et al. 
(2005), winds faster than predicted by CS06, or both. It would 
be worthwhile for other groups that possess hot Jupiter radia- 
tive transfer codes to compute radiative time constants at these 
temperatures and pressures. This is an area that we will pursue 
in the near future. 

Another possibility for a smaller planet-to-star flux ratio dur- 
ing secondary eclipse would be if the planet had a larger Bond 
albedo than calculations currently indicate. This would mean 
less absorbed stellar flux and correspondingly lower temper- 
atures everywhere on the planet. Cooler temperatures every- 
where on the planet would lead to chemical abundances that 
differ from our previous cases. More CH4 would form, at the 
expense of CO, which would also lead to a slightly higher mix- 
ing ratio for H2O, which shares oxygen with CO. 

Hot Jupiters are believed to have very low Bond albedos-on 
the order of 90% or more of incident stellar flux is expected to 
be absorbed. In Fortney et al. (2005), we found that our one- 
dimensional model atmosphere for HD 209458b scattered only 
8% of incident flux. For TrES-1, this was 6%. To date, there 
is at least a hint that the Bond albedo of TrES-1 may have been 
underestimated. Charbonneau et al. (2005), under the assump- 
tion that the planet emits as a blackbody, determined a Bond 
albedo of 0.31 ±0.14 from their Spitzer IRAC observations at 
4.5 and 8.0 pxn. Perhaps hot Jupiters are not quite as hot as had 
been previously thought. However, this determination should 
be considered very preliminary at this time. 

Richardson et al. (2006) have recently observed the transit 
of HD 209458b at 24 /im, as well. From an observed tran- 
sit depth of 0.0149 ± 0.0003, they determined the radius of the 
planet in this band to be 1.26 ±0.08 Rj, which includes uncer- 
tainties in the stellar radius. Our model predicts a change in 
the apparent planet-to-star flux ratio of ^0.00008 during the 20 
orbital degrees of the Richardson et al. (2006) transit observa- 
tions, ~4 times smaller than their uncertainty, and hence too 
small to have been detected. As was seen in Figures 9 and 10, 
the change in planetary flux as a function of orbital phase is not 
as pronounced near the time of transit (and secondary eclipse) 
as it is at other phases. 

5.3. Temperature Sensitivity 

To illustrate the sensitivity of our results to temperature 
changes, we have computed spectra as a function of orbital 
phase, using equilibrium (Case 0) chemistry for two additional 
models. These are additional dynamical models described in 
CS06, in which the base P-T profile of Iro et al. (2005) has 
been increased or decreased by 300 K, with correspondingly 
warmer or colder night sides. The full dynamical simulations 
have been run again with these parameters. The resulting light 
curves, in Spitzer bands, are shown in Figures 11a, and c. For 
the "cold" (-300 K) model CH4 is the dominant carbon car- 
rier on the night hemisphere, and CO on the day hemisphere, 
leading to large flux variation in excess of that found for our 
nominal CS06 simulation, as shown in Figure 1 1 . For the "hot" 
(+300 K) model, on both the day and night hemispheres, CO is 
the dominant carbon carrier. Methane absorption features are 
very weak on the night side, leading to flux variations in ev- 
ery band no larger than a factor of 2.7 from peak to trough. 
This model is somewhat similar to our earlier Case 2, but at 
warmer temperatures, as at all phases CO is the dominant car- 
bon carrier. Fluxes are everywhere greater in the hot model 
than the nominal model, and everywhere less in the cold model 



than in the nominal model. This is is due to the differing at- 
mospheric temperatures. The cold model best fits the Deming 
et al. (2005) datum at 24 ^m. In addition, the phase of maxi- 
mum and minimum flux in a given band can change by up to 
~10 orbital degrees between these simulations, due both to dif- 
fering chemistry and atmospheric dynamics. As discussed in 
Cooper & Showman (2005) and CS06, while these models all 
possess similarly strong east-to-west jets, the dynamical atmo- 
spheres differ slightly in detail. 

CS06 also examined non-equilibrium CH4/CO chemistry for 
these models. For the "cold" (-300 K) model, non-equilibrium 
chemistry was important, and the CH4/CO ratio at P < 1 bar be- 
came homogenized at 0.20 around the entire planet. This ratio 
was 0.014 for the nominal Case 1 decribed earlier. For the cold 
model equilibrium chemistry would predict a night side domi- 
nated by CH4 and a day side dominated by CO. For the "hot" 
model (+300 K) both equilibrium and non-equilibrium chem- 
istry predicts that CO is the dominant carbon carrier on both the 
day and side hemispheres. Figure lid shows our computed flux 
ratios for non-equilibrium chemistry for the cold case. As was 
shown previously for our nominal case, flux variation is smaller 
with non-equilibrium chemistry, because the mixing ratios of 
CH4, CO, and H2O are the same on the night and day hemi- 
spheres. Again, because of the nearly isothermal temperature 
structure of the day-side, chemical abundances have little effect 
on the day-side planet-to-star flux ratios, which are nearly equal 
for these two chemistry cases. 

These additional cases further highlight the importance of 
CH4/CO chemistry in the computation of infrared light curves. 
Infrared fluxes as a function of orbital phase are sensitive to 
the temperatures of the hemisphere facing the observer, as well 
as the abundances of important molecular absorbers. For plane- 
tary P-T profiles that cross important CH4/CO chemical bound- 
aries, as most hot Jupiters surely do, whether or not these 
species are found in their equilibrium mixing ratios has a major 
impact on resulting infrared flux in Spitzer bands, especially on 
the night side. 

5.4. The Future & Conclusions 

From the size of the error bar from the Deming et al. (2005) 
24 /im observations, it is clear that, could this instrument stabil- 
ity be sustained over the course of tens of hours of observations, 
the change in flux over time presented here could be detected. 
One-half of an orbital period for HD 209458b is 42 hours. If 
flux in the 8.0 /im band for HD 209458b is higher than models 
predict, as was the case for TrES-1 (Fortney et al. 2005; Bar- 
man et al. 2005), then this would be an attractive band as well, 
as the error bars should be smaller. 

What might one hope to see with sustained observations? If 
the peak in infrared flux does occur only ^25-30 orbital degrees 
before secondary eclipse, this would be difficulty to discern. 
However, the detection of any sort of ramp up in flux from the 
time of transit to secondary eclipse would give us important in- 
formation on the day-night temperature contrast. The recently 
discovered transiting planet HD 189733b (Bouchy et al. 2005) 
would likely be an even more attractive target, as the planet-star 
flux ratios are likely twice as large (Fortney et al. 2006; Dem- 
ing et al. 2006a), and the orbital period is ^40% shorter. We 
predict secondary eclipse planet-to-star flux ratios for this sys- 
tem in Fortney et al. (2006); our calculation at 16 /im is a good 
match to the published observation of Deming et al. (2006a). 

A clear prediction from our calculations here is that when 
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one uses realistic non-equilibrium chemistry calculations, the 
change in planetary flux as a function of orbital phase is greatly 
reduced, relative to equilibrium chemistry calculations, because 
the atmosphere's composition is fixed. For the Spitzer bands, 
for Case 0, the maximum flux variation is in the 3.6 /im band, 
which varies by a factor of 15 from peak to trough. The mini- 
mum variation is a factor of 2.2, in the 16 /im band. For Case 
1, this variation drops significantly, and the maximum variation 
is a factor of 3.7, in the 5.8 /im band, and the minimum is 2.0, 
again in the 16 /im band. We suggest that the 5.8 and 8.0 /im 
bands may be the best Spitzer bands in which to search for flux 
variations as a function of orbital phase, as these bands combine 
a high planet-to-star flux ratio and the sensitivity and stability 
of the IRAC detectors. 

If the day-side thermal emission of hot Jupiters is similar to a 
blackbody, as we find here, problems arise with the notion that 
the emission can be used to characterize the atmospheric chem- 
istry from secondary eclipse observations. Absorption features 
due to CH4, CO, and H2O would be nonexistent or extremely 
weak. Observations at other orbital phases would then take on 
additional importance. 

Given the significant spectral differences between our model 
and radiative-equilibrium models, it is clear that more work in 
this area is certainly warranted. We note that blackbody-like hot 
Jupiter emission can simultaneously explain all observations to 



date. First is the secondary eclipse observation at 24 /im by 
Deming et al. (2005), a clear detection. Second is the very 
low flux ratio upper limit in L' band by Deming et al. (2006b). 
Third is the 2.3 pm relative flux observation of Richardson 
et al. (2003a), which most one-dimensional models cannot fit 
(Seager et al. 2005). Fourth is the set of observations of TrES- 
1 by Charbonneau et al. (2005), who found an infrared spectral 
slope from 4.5 to 8 /im that was redder than that found by most 
one-dimensional models (Fortney et al. 2005; Barman et al. 
2005). Additional observations of these planets, especially in 
the Spitzer 3.6 /im band, which catches much of the predicted 4 
/im flux peak, would strengthen or refute this argument and pro- 
vide a critical test for the CS06 dynamical simulation. Deming 
et al. (2006b) posited, perhaps with a wink, that a blackbody- 
like spectrum would be more consistent with observations to 
date than any published hot Jupiter model. Due to the dynam- 
ically altered temperature structure of the atmosphere of HD 
209458b, we find that this could be the reality. 
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Table 1 
Chemistry Cases 



Case# 


XcH 4 


Xco 


Xh 2 o 


CH 4 /CO ratio 





equil. 


equil. 


equil. 


equil. 


1 


6.60E-6 


4.81E-4 


2.75E-4 


0.014 


2 


1.50E-8 


4.88E-4 


2.69E-4 


3.0E-5 



Note. — Mixing ratios are given for each chemistry 
case. Abundances from equilibrium chemistry calculations 
are used when "equil." is specified. 
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FIG. 1 . — Global temperature map at 2 mbar (top panel), 25 mbar (middle panel), and 220 mbar (bottom panel) for the CS06 dynamical simulation. Arrows show 
the direction and relative magnitudes of winds. Each longitude minor tick mark is 18 degrees and each latitude minor tick mark is 9 degrees. Each panel uses the 
same temperature shading scheme. 
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FIG. 2. — Top- view diagram of the the planetary orbit. The CS06 latitude/longitude grid covers the entire planet. The planet's night side is shaded gray. At every 
angle in the orbit (tp), the longitude is remapped to give the apparent longitude and latitude, as viewed by the distant observer. The cosine of the angle 9 made 
between a patch of the planet and the sub-observer point is /i. 
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FIG. 3. — In gray is the full pressure and temperature extent of the 3168 P-T profiles from the CS06 simulation. In thick black are 13 randomly selected profiles 
that allow one to make out individual profile shapes and temperatures. In the middle of the gray, shown in white long dashes, is the Iro et al. (2005) profile. The 
labeled dotted curves show the condensation curves of important cloud species. The dashed curve shows where CH 4 and CO have an equal abundance. All of these 
boundaries assume solar abundances. Note that day-side and limb profiles tend to be generally isothermal at P < 200 mbar. The night-side profiles tend to plot over 
one another on the far left. 
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FIG. 4. — The mixing ratios of CH4, CO, and H2O, as function of pressure, along three P-T profiles in the atmosphere of HD 209458b. These profiles are the 
one-dimensional profile from Fortney et al. (2005), in gray, the subsolar point of the CS06 simulation, in thin black, and the anti-stellar point of the CS06 simulation, 
in thick black. Arrows at the top of the plot indicate the fixed mixing ratios of Case 1 and 2 chemistry, given in Table 1. 
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FIG. 6. — Planetary emergent flux density (erg s _I cm -2 Hz -1 ) vs. wavelength as a function of orbital phase for Case equilibrium chemistry. The spectra are 
color-coded with the diagram in the upper right corner of the figure. All spectra are calculated at intervals of 60 orbital degrees. The black spectrum is for the night 
side of the planet, which is seen during transit. The red spectrum is 60 orbital degrees later. The dark blue spectrum shows the planet during secondary eclipse (180 
degrees), when the fully illuminated hemisphere is visible. The magenta spectrum is 300 orbital degrees after transit. The dashed black curve is the flux of a 1330 K 
blackbody, which plots behind the dark blue curve at A > 4/jm. Normalized Spitzer band passes are shown in dotted lines at the bottom and standard H, K, L, and 
M bands are shown at the top. 
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FIG. 7. — Planetary emergent flux density (erg s -1 cirT 2 Hz~') vs. wavelength at two orbital phases for all three of our chemistry cases. The spectrum at 180 
degrees, during the secondary eclipse, is essentially the same for all three cases. At zero degrees, during the transit, CH4 is very abundant in the Case trial, leading 
to deep absorption features centered at 1.6, 2.2, 3.3, and 7.8 fim. Absorption due to CO is muted in Case 0, but is strong in Cases 1 and 2. Water vapor absorption is 
strong in all cases. Inset: Zoom near 2.2 fim. The three thick horizontal lines connected by dashes are a constraint on the maximum relative height of the flux peak 
at 2.2 (im (middle band), compared to the average of the flux in the two adjacent bands, from Seager et al. (2005). The bars can be shifted together vertically. (See 
text.) 
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Fig. 8.— 



Effective temperature of the visible hemisphere of HD 209458b as a function of orbital phase for two chemistry cases. 
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FIG. 9. — Planetary flux in Spitzer bands as a function of planetary orbital phase. The left panels are Case and right panels Case 1. The bottom panels show the 
HD 209458b planet-to-star flux ratios. The top panels show the planetary flux at every phase divided by the planetary flux seen during the transit (at zero degrees). 
The dotted lines indicate predictions for Case 2 chemistry. At 180 orbital degrees is the Deming et al. (2005) secondary eclipse observation at 24 fjm, with one-sigma 
error bars (0.00260±0. 00046). This is shown in magenta, which is the color used for the 24 /im curve. 
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FIG. 10. — Planetary flux in MKO infrared bands as a function of planetary orbital phase. The left panels are Case and right panels Case 1. The bottom panels 
show the HD 209458b planet-to-star flux ratios. The top panels show the planetary flux at every phase divided by the planetary flux seen during the transit (at zero 
degrees). The dotted line indicates a prediction for Case 2 chemistry. At 180 orbital degrees is the Deming et al. (2006b) secondary eclipse observation at L' band, 
with one-sigma error bars. The observation is -0.0007±0.0014, so the actual data point is well off of the bottom of the plot. This is shown in dark blue, which is the 
color used for the L' curve. 
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FIG. 1 1. — Planetary flux in Spitzer bands as a function of planetary orbital phase. The lightcurves in panels A, B, and C assume equilibrium chemistry while 
panel D uses non-equilibrium chemistry. The bottom panels show the HD 209458b planet-to-star flux ratios. The top panels show the planetary flux at every phase 
divided by the planetary flux seen during the transit (at zero degrees). The left panel (A) shows the CS06 simulation with a 300 K increase in temperature. Panel 
B shows the nominal Case simulation, previously shown in Figure 9a. Panel C shows the CS06 dynamical simulation with a 300 K temperature decrease. The 
right-most panel (D) is also the simulation with the 300 K temperature decrease, but uses non-equilibrium chemistry with a fixed CH4/CO ratio of 0.20. (See text.) 
For all panels, the Deming et al. (2005) datum at 24 /xm is shown. 



